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We discuss the interior structure and composition of giant planets, and how this structure changes as 
these planets cool and contract over time. Here we define giant planets as those that have an observable 
hydrogen-helium envelope, which includes Jupiter-like planets, which are predominantly H/He gas, and 
Neptune-like planets which are predominantly composed of elements heavier than H/He. We describe 
the equations of state of planetary materials and the construction of static structural models and thermal 
evolution models. We apply these models to transiting planets close to their parent stars, as well as 
directly imaged planets far from their parent stars. Mechanisms that have been postulated to inflate the 
radii of close-in transiting planets are discussed. We also review knowledge gained from the study of the 
solar system's giant planets. The frontiers of giant planet physics are discussed with an eye towards future 
planetary discoveries. 



1. INTRODUCTION 

The vast majority of planetary mass in the solar system, 
and indeed the galaxy, is hidden from view in the interiors 
of giant planets. Beyond the simple accounting of mass, 
there are many reasons to understand these objects, which 
cut across several disciplines. Understanding the structure 
of these planets gives us our best evidence as to the forma- 
tion mode of giant planets, which tells us much about the 
planet formation process in general. As we will see, gi- 
ant planets are vast natural laboratories for simple materials 
under high pressure in regimes that are not yet accessible 
to experiment. With the recent rise in number and stunning 
diversity of giant planets, it is important to understand them 
as a class of astronomical objects. 

We would like to understand basic questions about the 
structure and composition of giant planets. Are they sim- 
ilar in composition to stars, predominantly hydrogen and 
helium with only a small mass fraction (~1%) of atoms 
more massive than helium? If giant planets are enhanced in 
"heavy elements" relative to stars, are the heavy elements 
predominantly mixed into the hydrogen-helium (H-He) en- 
velope, or are they mainly found in a central core? If a dense 
central core exists, how massive is it, what is its state (solid 
or liquid), and is it distinct or diluted into the above H-He 
envelope? Can we understand if a planet's heavy element 
mass fraction depends on that of its parent star? 

Giant planets are natural laboratories of hydrogen and 
helium in the megabar to gigabar pressure range, at temper- 
atures on the order of 10^ K. How do hydrogen and helium 
interact under these extreme conditions? Is the helium dis- 
tribution within a planet uniform and what does this tell us 



about how H and He mix at high pressure? What methods 
of energy transport are at work in the interiors of giant plan- 
ets? Can we explain planets' observable properties such as 
the luminosity and radius at a given age? 

The data that we use to shape our understanding of 
giant planets comes from a variety of sources. Labo- 
ratory data on the equation of state (EOS, the pressure- 
density-temperature relation) of hydrogen, helium, warm 
fluid "ices" such at water, ammonia, and methane, silicate 
rocks, and iron serve at the initial inputs into models. Im- 
portantly, data are only avaible over a small range of phase 
space, so that detailed theoretical EOS calculations are crit- 
ical to understanding the behavior of planetary materials at 
high pressure and temperature. Within the solar system, 
spacecraft data on planetary gravitational fields allows us 
to place constraints on the interior density distribution for 
Jupiter, Saturn, Uranus, and Neptune. For exoplanets, we 
often must make due with far simpler information, namely 
a planet's mass and radius only. For these distant planets, 
what we lack in detailed knowledge about particular plan- 
ets, we can make up for in number. 

Within six years of the Voyager 2 fly-by of Neptune, 
the encounter that completed our detailed census of the 
outer solar system, came the stu nning discoveries of the 
extrasolar giant planet 51 Peg b ( Mayor & Ouelod 1995 ) 
and also the first bo na fide brown dwarf, Gliese 229B 
(INakajima et alJll995h . We were not yet able to fully un- 
derstand the structure and evolution of the solar system's 
planets before we were given a vast array of new planets to 
understand. In particular the close-in orbit of 5 1 Peg b led 
to immediate questions regarding its history, structure, and 
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fate dGuillot et al.lll996l: iLin et al.lll996). Four years later 



the first transiting planet. HD 209458b (iCharbonneau et al 



tne nrst transiting planet^ i 

200(i iHenrv et alJ l200 Cn was found to have an inflated 
radius of ~1.3 Rj, confirming that proximity to a par- 
ent star can have d ramatic effects on planetary evolution 



dGuillot et al.lll99a) . The detections of over 50 additional 



transiting planets (as of August 2009) has conclusively 
shown that planets with masses greater than that of Sat- 
urn are composed predominantly of H/He, as expected. 
However, a great number of important questions have been 
raised. 

Much further from their parent stars, young luminous 
gas giant planets ar e being directly im a ged from the ground 
and from space dKalas et aL 2008 : Marois et al. 2008 



Lagrange et al. 20091) . For imaged planets, planetary ther 
mal emission is only detected in a few bands, and a planet's 
mass determination rests entirely on comparisons with 
thermal evolution models, that aim to predict a planet's 
luminosity and spectrum with time. However, the lumi- 
nosity of young planets is not yet confi dently understood 
dMarlev et al.ll20()7HChabrier et al.ll2007h . 

In this chapter we first outline the fundamental physics 
and equations for understanding the structure and thermal 
evolution of giant planets. We next describe the current 
state of knowledge of the solar system's giant planets. We 
then discuss current important issues in modeling exoplan- 
ets, and how models compare to observations of transiting 
planets, as well as directly imaged planets. We close with a 
look at the future science of extrasolar giant planets (EGPs). 



2. EQUATIONS AND MODEL BUILDING 

2.1 Properties Of Materials At High Pressure 

Materials in the interiors of giant planets are exposed 
to extreme temperature and pressure conditions reaching 
-10000 K and 100-1000 GPa (1-10 Mbai-). (See Figure 
1.) The characterization of materials' properties under these 
conditions has been one of the great challenges in experi- 
mental and theoretical high pressure physics. Ideally one 
would recreate such high pressure in the laboratory, char- 
acterize the state of matter, and then directly measure the 
equation of state (EOS) as well as the transport properties 
needed to model planetary interiors. While a number of 
key experiments have been performed, for a large part of 
Jupiter's interior we instead rely on theoretical methods. 

There are both static and dynamic methods to reach high 
pressures in laboratory experiments. The highest pressure 
in static compression experiments have been reached with 
diamond anvil cells ( Hemlev & Ashcrofti Il998t iHemlev 
20001: iLoubeyre et al.ll2002h . One has been able to reach 
'^400 GPa (4 Mbar), which exceeds the pressure in the cen- 
ter of the Earth but far less from —4000 GPa at Jupiter's 
core-envelope boundary. While many diamond anvil cell 
experiments were performed at room temperature, the com- 
bination with laser heating techniques has enabled one to 



approach some of temperature conditions that exist inside 
planets. 

The challenges in dynamic compression experiments 
are quite different. Reaching the required pressures is 
not the primary concern, but instead it is the difficulty in 
achieving a high enough density. Most dynamic exper- 
iments compress the material with a single shock wave. 
The locus of final p — T — P points reached by sin- 
gle shock from one particular initial point is called the 
Hugoniot . While this method pro vides direct access to 
the EOS ( Zeldovich & Raizei Il966h . the compression ratio 
rarely exceed values of 4. Instead the material is heated to 
very high temperature th a t exceed those in planetary inte- 
riors IjeanloretaDliQoS ImH^ Re- 
cently, static and dynamic compression techniques have 
been combined to address this issue (Eggert et al. 200l. 
By precompressing the sample in a diamond anvil cell be- 
fore a shock wave was launched, Eggert et al. were able to 
pr obe deeper into planetary interiors. Earlier experiments 
bv lWeir et aP dl996i) employed reverberating shock waves 
to reach high densities and thereby approached the state of 
metallic hydrogen at high temperature. 

The first laser shock experiments that reached megabar 
pressures predicted that the material to be highly com- 
pressible under shock conditions and to reach dens ities 
six times higher th an the initial state dda Silva et alj|l997l: 

— .,4. „i 1ir\r\oK tt — . i„+ :™ — Jt^ 1 .,4 



'Coninsetal"l998'). However, later experiments ("Knuds on et al 
2001; Knudson et al. 2003; Belov et al. 2002; Boriskov e tail 
|2005j) showed smaller compression ratios of about 4.3, 
whic h were in good agreeme nt with th eoretical predic- 
tions dLenoskvet al.|[r 997; Mili tzer & Cepe rlev 2000). 

Until very recently all models for giant planet interiors 
were based on chemical models ("Ebeling & Richert"l985'; 



Saumon & C habrier 1992; Saumon et al. 1995; Juranek & R edmei 
2000h that describe materials as a ensemble of stable 
molecules, atoms, ions, and free electrons. Approxima- 
tions are made to characterize their interactions. The free 
energy of the material is calculated with semi-analytical 
techniques and all other thermodynamic variables are de- 
rived from it. Chemical models require very little computer 
time, can easily cover orders of magnitude in pressure- 
temperature space, and have therefore been applied to nu- 
merous star and planet models. Chemical models do not 
attempt to characterize all the interactions in a many-body 
system. Molecular hydrogen at high density, e.g., is typi- 
cally approximated by a system of hard spheres where the 
excitation spectrum of the isolated molecule is modified by 
a density dependent term. In general, chemical models have 
difficulties predicting materials' properties in the strongly 
coupled regime where interaction effects dominated over 
kinetic effects. At high density and temperature where 
molecules are no longer stable or near a metal-to-insulator 
transition, chemical models require input from experiments 
or other theoretical techniques to fit adjustable parameters. 

Recently, progress in the field of theoretical description 
of planetary materials at high pressure has come from^rsf 
principles computer simulations. Such methods are based 
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Fig . 1 . — Density-temperature phase diagram of hot dense hydrogen. The dash-dotted hnes separate the molecular, atomic, 
metallic, and plasma regimes. The solid lines are isentropes for Jupiter and stars with 0.3, 1, and 15 solar masses. Single 
shock Hugoniot states as well as the inertial confinement fusion paths are indicated by dashed lines. The thin solid line 
shows p-T conditions of PIMC simulations. 



on the fundamental properties of electrons and nuclei and 
do not contain any parameters that are fit to experimental 
data. While approximations cannot be avoided altogether to 
efficiently derive a solution to the many-body Schrodinger 
equation, such approximations are not specific to a particu- 
lar material and have been tested for a wide range of mate- 
rials and different thermodynamic conditions. 

First-principles simulation can now routinely study the 
behavior of hundreds of particles at very different pressure 
and temperature conditions. Here we summarize three dif- 
ferent approaches: path integral Monte Carlo (PIMC), den- 
sity functional molecular dynamics (DFT-MD), as well as 
quantum Monte Carlo (QMC). The challenge of perform- 
ing accurate simulations has always been to make sure that 
the approximations, that are often necessary to perform the 
calculations at all, do not impact the predictions in a sig- 
nificant way. There are fundamental approximations such 
as the assumption of simplified functionals in DFT or the 
nodal app roximation in QMC and PIMC simulations with 
fermions (F oulkes et al.ll200lh . These approximations can 
in most cases only be checked by comparison between dif- 
ferent methods or with experimental results. Then there are 
also controlled approximations such as using a sufficiently 
large number of particles, long enough simulations, or a 
large enough basis set. These approximations can always be 
verified by investing additional computer time but it is not 
always possible to perform all tests at all thermodynamic 
conditions. 
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Fig. 2.— Comparison of the DFT-MD EOS with 
the Saumon-Chabrier (SC) and Sesame models. Three 
isotherms for pure hydrogen are shown in the metallic 
regime at high pressure. P is scaled by the volume per 
electron to the power | to remove most of the density de- 
pendence. 
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PIMC dPoUock & Ceperlevlll984l: ICeperlevlll995b is a 
finite-temperature quantum simulation method that exphc- 
itly constructs paths for electrons and nuclei. All correlation 
effects are included, which makes PIMC one of the most 
accurate finite-temperature quantum simulation methods 
available. The only fundamental approximation required is 
the fixed node approximation (Ceperley 1991, 1996) that is 
introduced to treat the fermion sign problem, which arises 
from the explicit treatment of electrons. The method is 
very efficient at high temperature and can provide one co- 
herent description of matter reaching up to a fully ion- 
ized plasma st ate. PI MC simulations have been applied 
to hydroge n (iPierleoni et al. 1994 ; Militzer & Ceperley 
200d |200lh. helium (lMihtzeili2006l l2009b . and their mix- 



occur gradually in the condit ion of giant planet interiors. 
Earlie r DFT-MD simulations ( Scandolol 2003 ; Bonev et al 



tures (' Militzer 2005V At low temperature, this method be 
comes more computationally demanding because the length 
of the path scales like l/T. At temperatures below ^5000 K 
where electronic excitations are not important, it is more ef- 
ficient to use a ground-state simulation method discussed 
next. 

Density functional molecular dynamics simulations rely 
on the Born-Oppenheimer approximation to separate the 
motion of electrons and nuclei. For a given configuration 
of nuclei, the instantaneous electron ground-state is derived 
from density functional theory. Forces are derived and nu- 
clei are propagated using classical molecular dynamics. Ex- 
cited electronic states can inc orporated by using the Mer- 
min functional ( Merminlll965 ). 

DFT is a mean field approach and approximations 
are made to treat electronic exchange and correlation ef- 
fects. Electronic excitations gaps are underestimated in 
many materials. A more sophisticated description of elec- 
tronic correlation effects is provided by quantum Monte 
Carlo jFoulkes et al. 2001 ) where one uses an ensemble of 
random walks to project out the ground-state wave func- 
tion. This method represents the ground-state analog of 
PIMC. To avoid the fermion sign problem, one also intro- 
duces a nodal approximation based on a trial wave func- 
tion. While most QMC calculations were performed for 
fixed nuclei, the method has recently been extended to 
fluids a nd calculations for fluid hydrogen have b een per- 
formed dPierleoni et alj|2004 lOelanev et al.ll2006l) . 

The EOS of dense hy drogen has been the subject of 



several DFT-MD studie s (ILenoskv et al.l 119971: iDesiarlais 



\2m5\ 'Bonev et al."2004'). Fig. |2] compares DFT-MD EOS 
from Militzer et al. (2008) with predictions from free en- 
ergy models. Even at the highest densities, one finds signif- 
icant deviations because no experimental data exist to guide 
free energy models. Furthermore, the density is not yet high 
enough for hydrogen to behave like an ideal Fermi gas. 

Of importance for the interiors of giant planets and 
the generation of their magnetic fields are the proper- 
ties of the insulator-to-me tal transition in dense hydro- 
gen dChabrier et al. l2006h . According to predictions 



from the best simulati on methods currentl y available. 



20041) had predicted a sharp dissociation transition but 
all these results have now been attributed to inaccuracies 
in the wave function propagati on with the Car-Parinello 
method dCar & Parrinelldll985l) . With the more accurate 



Born-Oppenheimer propagation meth od, the transition oc- 
curs g radually (see the discussion in Militzer & HubbardI 

('2009) ) but gives rise to a region of negative dP/ dT\\^ (" Bagnier et al 
2000; Vorbergeret al. 2007) that is shown in Fig. This 
leads to a negative Griineisen parameter and could intro- 
duce a barrier to convection in Jupiter and other giant plan- 
ets. However Militzer et al. (2008) demonstrated that in a 
hydrogen-helium mixture, the region of dP/dT\v < is 
shifted to lower temperatures than occur in Jupiter. 

The question whether the insulator-to-metal transition in 
hydrogen is smooth or of first order (plasma phase transi- 
tion) has been debated for a long time. A series of free 
ener gy models were constructed with a plasma phase transi- 
tion dEbeling et aPl 19931; ISaumon & Chabrieiiri992h. while 
others do not include one ("Ross" 1 998'; 'Sesame" 1992"). Re- 
cent shock wave experiments (Fortov et al. 2007) show evi- 
dence of an insulator-to-metal transition at temperature and 
pressure conditions that are consistent with theoretical pre- 
dictions. In their interpretation, the authors carefully sug- 
gest that their data may provide evidence for a first order 
phase transition. 



Rec ently Milit zer et al.l d2008h and independently lNettelmann et al 

d2008l) used first-principles simulations to derive an EOS 
to model Jupiter's interior While both studies relied on 
the same simulation technique, DFT-MD, the groups de- 
rived very different predictions for the size of Jupiter's 
rocky core, the distribution of heavy elements in its man- 
tle, and the temperature profile in its interior The dif- 
ferences between the two approaches are analyzed in 
Militzer & Hubbard ( 20091) and will be summarized later 
in this chapter. While some deviations in the simulation pa- 
rameters such as the number of particles and their effect on 
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Fig. 3. — Isochores derived from DFT-MD simulations of 
H-He mixtures (circles, F=0.2466) and pure hydrogen (tri- 
angles). Results for mixtures predict a positive Griineisen 
parameter, {dP/dT)\v > 0, along Jupiter's isentrope. 



MD dVorberger et al.l l2007h . this transition is expected to 
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the computed EOS have been identified, the differences in 
the predictions are mainly due to additional model assump- 
tions e.g. whether helium and heavy elements are homoge- 
neously distri buted throughout the in terior (.Militzer e t al., 
20081) . or not jNettelmann et al.ll2008l) . Only the deviations 
in the interior temperature profile is a direct consequence of 
the computed EOS. While pressure and internal energy can 
be obtained directly from simulations at constant volume 
and temperature, the entropy does not follow directly and is 
typicall y obtained by thermod yna mic integration of the free 
energy IMilitzer et alj (l2008b and iNettelmann etalj (|2008|) 
used different methods to compute the free energy but more 
work is needed to understand which method yields more 

accurate adiabats. 

Two recent papers jMorales et al. I I2OO9I: iLorenzen et al 



2009h provide evidence for the phase separation of the 
hydrogen-helium mixtures at the interior of Saturn and 
possibly also of Jupiter. Both papers rely on DFT-MD 
simulations but Morales et al. (2009) used larger and more 
accurate simulations and employed thermodynamic inte- 
gration to de termine the Gibbs free energy of mixing. 
Lorenzen et al. (2009) used a simplified approach where 
mixing entropy, which is the most difficult term to cal- 
culate, was taken from a noninteracting system of parti- 
cles. While others have^hown that non-ideal mixing ef- 
fects are important dVorberger et al.l l2007l) . the impact of 
this approximation on the pressure-temperature conditions 
where the hydrogen and helium start to phase-separate re- 
mains to be studied but it may explain whv lLorenzen et al.l 
("2009^ predict higher phase-separation temperatures than 
Morales et al. (2009). 

In addition to hydrogen and helium, water has also 



been extensively studied experi mentally (e.g. iLeeetal. 
2006 ) and computatio nally (e.g Cavazzoni et al.l Il999t 



Schwegler et all I2OO1I: Iprench et al.1 l2009l) . Often water, 
ammonia, and methane, are grouped together as "plane- 
tary ices," as a generic phrase for 0-, C-, and N-dominated 
volatiles (mostly H2O, NH3, and CH4), which are likely 
found in fluid form, not solid, within giant planets. Fur- 
thermore, these components are probably not found as in- 
tact molecules at high pressure. Water is found from first- 
principles calculations to dissociate into HgO ^ + OH ~ ion 
pairs above ~2000 K at 0.3 Mbar dC avazzo ni et al. 



Schwegler elal] I2OOII: iMattsson & Desiarlaisrf2006l) . and 



19991 



indeed high electr ical conductivitie s are measured at pres- 
sures near 1 Mbar JChau et alJ200ll) . Moving to even heav- 
ier elements, "rock" refers primarily to silicates (Mg-, Si- 
and O-rich compoun ds) and often includes iron and other 
"metals" as weU (see IStevensonll 1 9851; iHubbardI 1 984 . Un- 
certainties in the EOSs of heavy elements are generally less 
important in structure and evolution calculations than those 
for hydrogen and heliu m, but they certainly do have impor- 
tant q uantitative effects ( Hubbard et al Jl 1 99l : Baraffe et al 
20081) . 



This section is only able to give a flavor of the vast 
array of EOS science going on at the boundary between 
physics and planetary sciences. We are now in an era 



where long sought-after advances in experiment and first- 
principles theory are occurring at a fast pace. 

2.2 Basic Equations 

To a first approximation, a planet can be considered as a 
non-rotating, non-magnetic, fluid object where gravity and 
gas pressur^H are the two main contributors to forces in 
the interior. Planetary structure and evolution can thus be 
described in a spherically symmetric configuration and are 
governed by the following conservation equations: 



Mass conservation 



dm 2 
— = Airr p, 



Hydrostatic equilibrium 



dP 
dr 



(1) 



(2) 



dL dS 
Energy conservation : -— = —iTrr^pT-—, (3) 
or at 

where L is the intrinsic luminosity, i.e the net rate of radial 
energy flowing through a sphere of radius r. The rate of 
change of the matter entropy is due to the variation of its 
internal energy and to compression or expansion work, ac- 
cording to the first and second laws of thermodynamics. A 
complete model requires an additional equation describing 
the transport of energy in the planet: 



dT 
dr 



dPT 
'd^P ' 



(4) 



where V = is the temperature gradient. If energy 

transport is due to radiation or conduction, it is well de- 
scribed by a diffusion process with: 



dT 

dr 



kLP 



IdnacG mT'^ 



(5) 



The total opacity of matter, k, accounts for radiative and 
conductive transport and is defined by — Kross + '^cond 
with kross and Kcond are the Rosseland mean radiative 
opacity and the conductive opacity, respectively. Energy 
transport by conduction in a fluid results from collisions 
during random motion of particles. For planets essentially 
composed of H/He, energy transfer is due to electrons in 
the central ionized part, whereas molecular motion domi- 
nates in the outer envelope. Because of the high opacity of 
H/He matter in planetary interiors, convection is thought to 
be the main energy transport and the temperature gradient 
is given by the adiabatic gradient Vad — ( ^;"p )s ■ If the 
planet has a core made of heavy material (water or ice, rock, 
iron), heat transport can be due to convection or conduction 
(electrons or phonons), depending on the core material, its 
state (solid or liquid) and the age of the planet. 

2.3 Thermal Evolution And Atmospheric Boundary 
Conditions 



^ Radiation pressure is negligible in planetary interiors. 
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Starting from a high entropy, hot initial state, the lumi- 
nosity of a planet during its entire evolution is powered by 
the release of its gravitational Eg and internal Ei energy and 
is given by (see Eq. (3)): 



L{t) = -j^{Eg+E,) = 



„ d , r , f de ^ 
P—( — jdm — I —dm, 



M 



dt p 



M 



dt 



(6) 

where e is the specific internal energy. The virial theorem, 
which applies to a self-gravitating gas sphere in hydrostatic 
equilibrium, relates the thermal energy of a planet (or star) 
to its gravitational energy as follows: 



aEi + = 0, 



(7) 



with a=2 for a monoatomic ideal gas or a fully non- 
relativistic degenerate gas, and a = f for an ideal diatomic 
gas. Contributions arising from interactions between parti - 
cles yield corrections to the ideal EOS (see lGuillotil2005b . 
The case a = | applies to the molecular hydrogen outer 
regions of a giant planet. Note that the mass fraction in- 
volved in these regions, for a Jupiter-like planet, is usually 
negligible compared to that involved in the central core and 
the metallic H region. 

According to Eqs. (6) and (7), the planet radiates: 



L oc 



dEg 
dt 



defining a characteristic thermal timescale tkh^ 



E. 

T 



GAP 
RL 



(8) 



(9) 



For a 1 Mj gaseous planet, with negligible heavy element 
content, tkh ~ 10^ yr at the beginning of its evolution r| 
and TKH > 10^° yr after 1 Gyr (see e.g. lBaraffe et al.l2003l 
for values of R and L at a given age). The reader must keep 
in mind that Equation (9) is a rough estimate of the char- 
acteristic timescale for cooling and contraction of a planet. 
This timescale can be longer by 1 or 2 orders of magni- 
tude than the value derived from Eq. (9). Indeed, when 
degeneracy sets in, a significant fraction of the gravitational 
energy due to contraction is used to increase the pressure 
of the (partially) degenerate electrons and the luminosity of 
the planet is essentially pr ovided by the thermal cooling of 
the ions (see lGuillotll2005b . 

The rate at which internal heat escapes from a planet 
depends on its atmospheric surface properties, and thus 
on the outer boundary conditions connecting inner and at- 
mospheric structures. Put another way, although the inte- 
rior may be efficiently convecting, the radiative atmosphere 
serves as the bottleneck for planetary cooling. For objects 
with cold molecular atmospheres, the traditional Edding- 
ton approximation assuming that the effective temperature 
equals the local temperature at an optical depth tross ~ 



^Note that this value is highly uncertain since it depends on the initial state 
and thus on the details of the planet fomiation process. See §4.6 



2/3 pro vides incorrect therrnal pro files and large errors on 
7c ff (see IChabrier & Baraffell2000l and references therein). 
Modern models for planets incorporate more realistic atmo- 
spheric boundary conditions using frequency dependent at- 
mosphere codes. Inner and outer temperature-pressure pro- 
files must be connected at depths where the atmosphere be- 
comes fully convective, implying an adiabatic thermal pro- 
file, and optical depth is greater than one. Th e connection is 



done at a fixed pressure, usually a few bars ( Burro ws et al 



[1997, 2003; Fortney & Hubbard 2003; Guillot 2005i) or at a 
fixed optical depth, usually at r = 100 dChabrier & Baraffe 
I2OO O). The numerical radius corresponding to the outer 
boundary conditions provides, to an excellent approxima- 
tion, the planet's photospheric radius, where the bulk of the 
flux escapes. 

2.4 Effect Of Rotation 

The solar system giant planets are relatively fast rota- 
tors with periods of about 10 hours for Jupiter and Saturn, 
and a bout 17 hours for Neptune and Uranus (see iGuillot 



20051 and references therein). Rotation modifies the inter- 



nal structure of a fluid body and yields departures from a 
spherically symmetric configuration. Its effects can be ac- 
counted for using a perturbation theory, which has been 
extensively de veloped during the past century for plan- 



ets and stars (iTassoull 119781 ; IZharkov & Trubitsvnl 11974 , 
Il976. 1978 ). An abundant Uterature exists on the appli- 
cation of th is theory to models for our solar system gi- 
ant planets dHubbard & Marleylll989l iGuillot et al.1ll994 : 
[Podolak et al.lll995[) . One often refers to the so-called the- 
ory of figures, presented in details in the works of Zharkov 
and Trubitsyn. 

Here we only briefly explain the concept of the the- 
ory. The underlying idea is to define a rotational poten- 
tial V^iot such that surfaces of constant P, p and U coin- 
cide. U = Vrot + Vgrav IS the total potential (rotational -i- 
gravitational) and obeys the hydrostatic equilibrium equa- 
tion VP = pVU . The solution for the level surfaces (or 
figures) can be expressed in terms of an expansion in even 
Legendre polynomials, also called zonal harmonics, P2n- 
As a consequence of rotation, the gravitational potential de- 
parts from a spherically symmetric potential and can be ex- 
pressed in terms of even Legendre polynomials and grav- 
itational moments J2n (also called zonal gravitational har- 
monics) as follows: 



GM 



OO 

l-E(7)'""^2nP2n(cOS0) 
n— 1 



(10) 



where r is the distance to the planet center, a its equatorial 
radius and the polar angle to the rotation axis. The gravi- 
tational moments can be related to the density profile of the 
planet (see Zharkov & Trubitsyn, 1974) as follows: 



J2i\ 



1 



p{r,6yP2,,{cos6)dv, (11) 
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where dv is a volume element. A measure of J2n by 
a spacecraft coming close to a planet thus provides con- 
straints on its density profile. This powerful and elegant 
method was applied to our four giant planets using the de- 
termination of J2, J4 and Jg from the trajectories of the 
space missions Pioneer and Voyager (see Guillot, 2005, for 
a review, and the section below, "Solar System Giant Plan- 
ets"). 

As is well known from the theory of figures, the low or- 
der gravitational moments provide constraints on the den- 
sity/pressure profile in the metallic/molecular hyd rogen re- 
gions but do not sound the most central regions (IHubbard 
I999I see also Fig. 4 of Guillot, 2005). Consequently, the 
presence of a central core, its mass, composition and struc- 
ture, can only be indirectly inferred from the constraints on 
the envelope provided by the gravitational moments. 

2.5 Mass-Radius Relation 

The mass-radius relationship for planets, and more gen- 
erally for substellar/stellar objects, contains essential in- 
formations about their main composition and the state 
of matter in their interior. The fundamental work by 



electron degeneracy effects and renders the radius almost 
constant with mass around the critical mass. The full cal- 
culation yields, for gaseous H/He planets, Afcrit ~ 3Afj, 
amazingl y close to the results based o n the simplified ap- 
proach of IZapolskv & Salpeten i 



19691) . Below the critical 



Zapolsky & Salpeterl (Il969h is a perfect illustration of this 
statement. The analysis of cold (zero-temperature) spheri- 
cal bodies of a given chemical composition and in hydro- 
static equilibrium shows the existence of a unique mass- 
radius relation and of a maximum radius i?niax at a critical 
mass Merit- The very existence of a maximum radius stems 
from two competing physical effects characteristic of the 
state of matter under planetary conditions. The first effect 
is due to electron degeneracy, which dominates at large 
masses and yields a mass-radius relationsh ip R ix M~^^^ 
chara cteristic of fully degenerate bodies (IChandrasekhar 
1939h . The second effect stems from the classical elec- 
trostatic contribution from ions (Coulomb effects) which 
yields a mass radius relation R oc M^/^, characteristic 
of inco mpressible Earth-like planets. IZapolskv & Salpetei 
(Il969h ind a critical mass of 2.6 Afj where the radius 
reaches a maximum value i?max ~1 ^,1 for a gaseous H/He 
planet. The critical mass increases as the heavy element 
content increases, while i?inax decreases. 

The true mass-radius relationship, derived from models 
taking into account a realistic equation of state (see §2.1) 
yields a smoother dependence of radius with mass, as dis- 
played in Fig. |4] The transition between stars and brown 
dwarfs marks the onset of electron degeneracy, which in- 
hibits the stabihzing generation of nuclear energy by hy- 
drogen burning The t ypical transition mass is ~ 0.07 Mq 
dBurrows et all I2OOII: IChabrier & Baraffd I2OOOI) . Above 
this transition mass, the nearly classical ideal gas yields 
a mass-radius relationship R oc M. In the brown dwarf 
regime the dominant contribution of partially degenerate 
electrons, balanced by the contribution from ion interac- 
tions yields R oc AI^^^^ instead of the steeper relation- 
ship for fully degenerate objects. The increasing contribu- 
tion of Coulomb effects as mass decreases competes with 



mass. Coulomb effects slightly dominates over partially de- 
generate effects, yielding a smooth variation of radius with 
mass close to the relation R oc Af 



3. SOLAR SYSTEM GIANT PLANETS 



3.1 Jupiter And Saturn 

Jupiter and Saturn, composed primarily of hydrogen and 
helium, serve as the benchmark planets for our understand- 
ing of gas giants. Indeed, these planets serve as calibrators 
for the structure and cooling theory used for all giant plan- 
ets and brown dwarfs. The significant strides that have been 
made in understanding these planets have come from theo- 
retical and experimental work on the equation of state of 
hydrogen, as well at the spacecraft-measured gravitational 
moments. 

As discussed in S|2] with a basic knowledge of the EOS 
of hydrogen, helium, and heavier elements, along with ob- 
servations of Jupiter's mass and radius, one can deduce 
that Jupit er and S aturn are mostly composed of hydrogen 
(see, e.g. [Seager et al. 2007). The same argument shows 
that Uranus and Neptune are mostly composed of elements 
much heavier than hydrogen. 

Models of the interior structure of Jupiter in particular 
were investigated by many authors in the 20th century, and 
investigations of Jupiter and Saturn became more frequent 
with the rise of modern planetary science in the 1950s and 
1960s. However, even these initial models were uncertain 
as to whether the hydrogen is solid or fluid. The path to- 
wards our current understan ding o f giant planets started 
with the observation bv lLow ( 1966h that Jupiter emits more 
mid infrared flux tha n it received from the Sun. Soon af- 
ter, iHubbardI d 1968b showed that the observed heat flux 



could only be carried by convection (as opposed in radia- 
tion or conduction), showing that Jupiter's interior is fluid, 
not solid. Thus began the paradigm of giant planets as hy- 
drogen dominated, warm, fluid, convective objects. 

However, it was clear that Jupiter and Saturn were 
not composed entirely of hydrogen and helium in so- 
lar proportions — the y have radii too small for their mass. 
Podolak & Cameroiil ( 1974 ) showed that these planets 
likely had massive heavy elements cores, which tied 

their interior structure to a possible format i on mechanism 

JPerri & Cameronll974l:lMizuno et alJl978l) . [stevenson & Salpetei 
(1977a^ performed detailed investigations of phase dia- 
grams and transport properties of fluid hydrogen and he- 
lium, putting our understanding of these planets on a firmer 
theoretical footing. The issues discussed in these papers 
still reward detailed study. Also in the early to mid 1970s, 
the first thermal evolution models of Jupiter were computed 
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jGrossman et al. 1972 ; Bodenheimeill 1 974 ; Graboske et al . 
1975h and connections were made between these models 



and similar kinds of cooling models for fully convective 
very low mass stars. Much of the early kn owledge of gian t 
pla net structure and evolution is found in iHubbardI (119731) 
and Zharkov & Trubitsvnl(ll978b . A consistent finding over 
the past few decades is that fully adiabatic cooling models 
of Jupiter can reproduce its current Tcff to within a few K. 

The situation for Saturn is more complicated. There is 
a long-standing cooling shortfall in thermal evolution mod- 
els of Saturn: the planet is ~50% more luminous than one 
calculates for a 4.5 Gyr-old, adiabatic, well-mixed planet 
JPollack et all Il977t [stevenson & Salpeter 1977a). The 
most likely explanation is that the He is curre ntly phase 
separ ating from the liquid metallic hydrogen ([Stevenson 
19751) . and has been for the past 2-2.5 Gyr. This immis- 



cible He should coalesce to form droplets, that are denser 
than the surrounding H/He mixture, and then "rain" down 
within the planet. This differentiation is a change of gravita- 
tional potential energy into thermal energy. Recent models 
of the evolution of Saturn with the additional energy source 
indic ate that the He may be r aining down on top of the 
core (iFortney & Hubbardll2003h . The largest uncertainly in 
properly including He phase separation into cooling mod- 
els is the current observed abundance of He is Saturn's at- 
mosphere, which does appear to be depleted in He rela- 
tive to pr otosolar abundances, but the error bars are large 
JConrath & Gautier 2000)0 The in situ observation that 
Jupiter's atmosphere is modestly depleted in helium as well 
shows that our understanding of the evolution of both these 
planets is not yet complete. 

Figure [5] sh ows models of Jupiter's and Saturn's interi- 
ors taken from lGuilloi (2005). As Jupiter is 3.3 times more 
massive than Saturn, a greater fraction of its interior mass 
is found at high pressure. Therefore, most of the hydrogen 
mass of Jupiter is fluid metallic, while for Saturn it is fluid 
molecular, H2. Jupiter's visible atmosphere is warmer than 
Saturn's. If their interiors are fully adiabatic, this implies 
that Jupiter is always warmer at a given pressure than Sat- 
urn. Since Jupiter has a larger interior heat content (residual 
energy leftover from formation), this also means that at a 
given age, Jupiter is always more luminous than Saturn. 

Structural models for Saturn show that it is more cen- 
trally condensed than Jupiter Although significant cen- 
tral condensation is expected due to the compressibility of 
H/He, detailed models of Saturn indicate even greater cen- 
tral condensation is needed to match its gravitational field. 
This confirms that a significant fraction of the heavy ele- 
ments must be in the form of central de nse core. Recent es- 
timates from Saumon & GuillotI ( 20041) indicate a core mass 
of 10-20 M(^. The majority of Saturn's heavy elements are 
within the core. For Jupiter, the situation is less clear cut for 
several reasons. First, a 10 core would only be 3% of 



^Note that this "He rain" cannot power the inflated hot Jupiters, as their 
large radii require interior temperatures warmer than needed for phase sep- 
aration. 
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Fig. 5. — I nterior views of Jup iter and Saturn, from calcula- 
tions with the lSaumon et alj ( ll995h "chemical picture" H/He EOS. 
Jupiter is more massive, which leads to a greater fraction of its 
mass in the high-pressure liquid metallic phase. It also has a higher 
temp erature at a gi ven pressure. Interior temperatures are taken 
from lGuillotlbOOSi) . 
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the planet's total mass, so EOS inputs must be accurate to 
this same percentage for real constraints on the core mass. 
Second, since Jupiter is more massive than Saturn, a greater 
fraction of its inter ior is in the 1- 100 Mbar region that it 
difficu lt to model. iGuillotI (Il999h and ISaumon & GuillotI 
(120041) found that Jupiter models without any core were al- 
lowed, and that the majority of the planet's heavy elements 
must be mixed into the H/He envelope. This would be a 
clear difference compared to Saturn. On the whole Sat- 
urn possesses < 25 of heavy elements, while Jupiter 
is < 40 Mq. Saturn appears to be a relatively greater frac- 
tion of heavy elements by mass. 

Very recently, two groups have computed new models of 
the interior of Jupiter, based on firs t-principles equations 
of state for hydrogen and helium. iMilitzer et al.l ( 2008 ) 



simulate the hydrogen-helium mixtures directly, under the 
pressure-temperature conditions found in the interior of 
Jupiter. This work predicts a large core of 14 - 18 for 
Jupiter, which is in line with estimates for Saturn and sug- 
gests that both planets may have formed by core-accretion. 
The paper further predicts a small fraction of planetary ices 
in Jupiter's envelope, suggesting that the ices were incor- 
porated into the core during formation rather than accreted 
along with the gas envelope. Jupiter is predicted to have 
an isentropic and fully convective envelope that is of con- 
stant chemical composition. In order to match the observed 
gravitational moment J4, the authors suggest that Jupiter 
may not rotate as a solid body and predicted the existence 
of deep winds in the interior leading to differential rota- 
tion on cylinders. Differential rotation as well as the size of 
Jupiter's core can potentially be measured by the forthcom- 
ing Juno orbiter mission, briefly described below. 

Alternativelv. i Nettelmann et all ( I2OO8I) computed EOSs 
for hydrogen, helium, and water separately, and investi- 
gate interior models under the assumption of an additive 
volume rule at constant temperature and pressure. They 
hypothesize a heavy element-enriched metallic region and 
heavy element-depleted molecular region. They overall 
find a large abundance of heavy elements within Jupiter, 
but far less in a distinct core. A further key difference 
between these two models that has not yet been investi- 
gated is interior temperatures-the INettelmann et al. I (l2008h 



model predicts higher temperatures than do iMihtzer et al 



(120081) . This will have important consequences for cooling 
mo dels of the planet, which will need to be investigated. 
The INettelmann et all (l2008h and flie M iUtzer et al. (200l) 
approaches are compared in iMilitzer & HubbardI (l2009b . 
while Figure |6]gives a schematic comparison. 

The resolution of the H/He phase separation issue in Sat- 
urn may only come from an entry probe into Saturn's atmo- 
sphere, which could measure the He/H2 ratio in situ, as was 
done for Jupiter. If the He mixing ratios were known pre- 
cisely in each atmosphere, it is possible that phase diagram 
and evolution history combination could be derived that 
constrains the demixing region of the hig h-pressure H/He 
phase diagram (iFortney & HubbardlEool . In theory, if at- 
mospheric abundances of He, C, N, and O were known, this 



would be a strong constraint on the properties of the en- 
tire H/He envelope. In practice, this has not yet happened 
for Jupiter because the Galileo Entry Probe was apparently 
only able to measure a lower limit on the O abundance from 
water vapor The forthcoming Juno mission will use mi- 
crowave spectroscopy to measure the deep H2O and NH3 
abundances at pressures near 100 bars (below cloud con- 
densation levels), which should help to resolve this issue 
for Jupiter Juno is expected to launch in 2011 and arrive 
at Jupiter in 2016, with detailed data analysis beginning in 
2018. 

Juno, which is a low periapse orbiter, will also exquisitely 
map the planet's gravitational and magnetic fields. Work by 
HubbardI d 1999b has shown that if surface zonal flows ex- 



tend down to 1000 km depth (P ^ 10 kbar) then this should 
be observable in the planet's gravity field. This will give us 
a view of the mechanics of the interior of the planet. Does 
it rotate as a solid body, on cylinders, or something more 
exotic? Furthermore, the tidal response of Jupiter to the 
Galilean satellites may also be detectable, which will help 
to constrain the planet's core mass. 

The Juno orbiter will also map Jupiter's magnetic field. 
The magnetic fields of Jupiter and Saturn are large and pre- 
dominantly dipolar, with a small tilt from the planet's rota- 
tion axis, similar to the Earth. This tilt is 9.6° for Jupiter, 
but less then 1° for Saturn. These fields are consistent with 
their production via a dynamo mechanism within the liq- 
uid metallic region of the interior Unfortunately, dynamo 
physics is not understood well enough to place strong con- 
straints on the interiors of Jupiter and Saturn. 

3.2 Uranus And Neptune 

Uranus and Neptune have not received the same atten- 
tion that have been paid to Jupiter and Saturn. The neglect 
of the "ice giants" has been due to relatively less precise ob- 
servational data, and the complicated picture that the data 
has revealed. The first stumbling block has been the rel- 
atively high density of these planets, compared to Jupiter 
and Saturn. This high density shows that these planets are 
not predominantly composed of hydrogen and helium. But 
then what is their composition? Mostly the fluid planetary 
ices of water, ammonia, and methane? Mostly rock and 
iron? The mass/radius of these planets can be matched 
with a very wide range of compositions f rom the three 



categories of H/He gas, ice, and rock (e.g. iPodolak et al 
19951). A high pressure mixture of rock and H/He can 



very nicely mimic th e pressure-density relation of the ices 
dPodolak et al.lll99lh . which means that even gravity field 
data, which helps to elucidate central condensation, cannot 
break the ice vs. rock/gas degeneracy. Therefore, modelers 
often have had to resort to cosmogonical constraints, such 
as an assumed ice-to-rock ratio from protosolar abundance 
arguments, to constrain structural models. Recent estimates 
of the H/He mass fraction yield 1-2 in both planets 
(with an ice-to-rock ratio of 2.5), with a hard upper limit 
of in each if only rock and H/He gas are present 
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jPodolak et alJbOOOl) . 

Nearly all of collected knowledge on the structure and 
evolution of these planets have been gathered in the Uranus 
and Neptun e and Triton Arizona Spa ce Science Series boo k 
chapters bv lPodolak erall (Il99lb and lHubbard et alj (1 1995b . 
respectively. These chapters summarize our state of knowl- 
edge of Uranus and Neptune as of the early 1990s, after the 
Voyager 2 encounters. Importantly, novel research is also 
presented in those chapters that is not found in the more 
readily available literature. Both chapters are worth de- 
tailed study. Since that time, the only more recent update is 
the work of Marley and Podolak, who investigated Monte 
Carlo interior models of these planets with a minimum of 



assumptions regard ing interior density (IMarley et al.lll995 



Podolak et al.l 1200 0). Without assumptions regarding the 



layering of gas, ices, and rock, these authors perform Monte 
Carlo studies of the interior density distribution, which uses 
the gravity field alone to show where density jumps, if any, 
must occur 

The second major complication with these planets, after 
composition degeneracy, is the interior heat flow. While 
fully adiabatic, fully convective thermal evolution mod- 
els reproduce the current luminosity of Jupiter, and un- 
derpredict the luminosity of Saturn, they overpredict the 
luminosity of Neptune and Uranus. The situation for 
Uranus is especially dramatic, as no intrinsic flux from 
the planet's interior was detected by Voyager 2. At least 
two important ideas partially address the heat flow issue. 
Hubbard & Macfarland (fl980h suggested that the absorbed 
and reradiated stellar flux may be large enough to swamp 
the intrinsic flux, a smaller component. This would be 
a larger effect in Uranus than Neptune, since it is closer 
to the Sun. This same effect, on a much more dramatic 
scale, is seen for the hot Jupiters, where the intrinsic flux 
is unmeasurable, since it is 10^ smaller than reradiated ab- 
sorbed flux. While this effect is certainly real in Uranus 
and Neptune, it alon e cannot explain the low heat flows 
(lHubbardetal.lll995h . 

The problem may well be in the assumption that the in- 
terior is partitioned into well defined layers of H/He gas, 
the fluid ices, and rock. If these distinct layers exist, then 
convection should be efficient in each layer, and the interior 
heat should be readily transported to the surface. However, 
it is well known that composition gradients can readily sup- 
press convection, as a much steeper temperature gradient is 
needed for convective instability to occur, from the Ledoux 
criterion. If large regions of the interior of these planets are 
stably stratified, then stored residual energy from forma- 
tion will be "locked" into the deep interior, and will only be 
transported quite slowly. At gigayear ages this would lead 
to a small intrinsic luminosity. A promising explanation for 
the reduced heat flow of Uranus and Neptune is that the 
deep interiors of the planets, which are likely a mix of fluid 
ices and solid rock, are predominantly stratified, with only 
the outer ~ 1/3 of the heavy elem ent interior region freely 
convecting (iHubbard et al ] [l995l) . 

Recently, interior geometry as proposed above was in- 
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stable due to composition gradients. Circles with arrow heads in- 
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fraction of heavy elements and it may have a larger freely convec- 
tive region. This and other inferences are, however, uncertain. 
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vestig ated with 3D numerical dynamo models dStanlev & Bloxhhflij rradiation Effects 



2006h . Uranus and Neptune are known to have complex 
magnetic fields that are non-dipolar, non-axisymmetric, and 
tilted significantly from their rotation axes. Stanley & Blox- 
ham have found that they can reproduce the major features 
of these fields with an outer convecting ionic shell, in the 
outer '^20-40% of the interior region, with the innermost 
regions not contributing to the dynamo. This is strong ev- 
idence that the interiors of these planets are complex and 
are not fully convective. In addition, structure models that 
include a pure H/He envelope, a pure icy layer, and an in- 
ner rocky core are not consistent with the gravity field data 
of either planet. A possible interior view of each planet is 
shown in Figure [T] 

Since we will not have additional data to constrain the 
gravitational moments of these planets, or an entry probe, 
for perhaps decades, further progress in understanding the 
structure and evolution of Uranus and Neptune must come 
from new theoretical ideas. At this time, new thermal evo- 
lution models are needed to explore in some detail what 
regions of the interior are indeed convective. Future work 
on heat transport in double-diffusive convective regions in 
planetary interiors (see |j4]l should focus on Uranus and 
Neptune as well as exoplanets. 



4. RECENT HIGHLIGHTS: TRANSITING AND 
DIRECTLY IMAGED PLANETS 



4.1 The Observed Mass-Radius Relationship 

The determination of mass-radius relationships of exo- 
planets, using photometric transit and Doppler follow-up 
techniques, provides an unprecedented opportunity to ex- 
tend our knowledge on planetary structure and composition. 
The 20th century was marked by the solar system explo- 
ration by space missions, revealing the complexity and di- 
versity of giant planets in terms of chemical composition 
and internal structure. The prolific fishing for transiting ex- 
oplanets in the beginning of this 21st century, with a catch 
of over fifty objects, not only confirms this diversity but 
also raises new questions in the field of planetary science. 
Two benchmark discoveries illustrate the surprises planet 
hunters were faced with. The very first transiting planet 
ever dis covered, HD 209458b dCharbonneau et all l2000t 
Henrv et al. 2000), was found with an abnormally large ra- 
dius, a puzzling property now shared by a growing frac- 
tion of transiting exoplane ts. At the other e xtreme, a Saturn 
mass planet, HD 149026b dSato et al.ll2005b was discovered 
with such a small radius that more than 70 M® of heavy 
elements is required to explain its compact structure. This 
discovery raised in particular new questions on the forma- 
tion process of planets with such a large amount of heavy 
material. The diversity in mean density of transiting planets 
yet discovered is illustrated in Fig. [8] 



The discovery of HD 209458b and of the additional tran- 
siting exoplanets that followed has opened a new era in gi- 
ant planet modeling. The modern theory of exoplanet radii 
starts with models including irradiation effects from the par- 
ent star These effects on planet evolution are accounted for 
through the coupling between inner structure models and 
irradiated atmosphere models, following the same method 
described in §2.3. Current treatments are based on simpli- 
fied treatments of the atmosphere, using ID plane-parallel 
atmosphere codes. They however allow one to understand 
the main effects on planetary evolution. In reality the im- 
pinging stellar flux has an angle of incidence which is a 
function of the latitude and longitude, but in ID one at- 
tempts to compute a planet-wide or day-side average at- 
mosphere profile, using a parameter / which represents the 
redistribution factor of the stellar flux over the planet sur- 



face (Baraffe et al. 2003; Burrows et al. 2003; Fortn ev et al 



l2006h . The incident stellar flux Fine is explicitly included 
in the solution of the radiative transfer equation and in the 
computation of the atmospheric structure and is defined by: 



\nc. 



(12) 



where i?» and are the stellar radius and flux respec- 
tively, and a the orbital separation. The current generation 
of models often use /=1, corresponding to a stellar flux re- 
distributed over the entire planet's surface or / = 2 if heat 
is redistributed only over the day side. Heat redistribution 
is a complex problem of atmospheric dynamics, depend- 
ing in particular on the efficiency of winds to redistribute 
energy from the day side to the night side. This question 
is a challenge for atmospheric circulation modelers (Show 



man et al. chapter, this volume, as well as I Showman et al 



2008bl) . This nascent field is growing rapidly with observa- 



tional constraints provided by infrared lightcurves obtained 
with Spitzer, which are starting to provide information on 
the temperature structure, composition and dynamics of ex- 
oplanet atmospheres. 

Although treatments of irradiation effects can differ in 
the details, with possible refinements accoun ting for phase 
and a n gle dependences o f the in cident flux (Barm an et alJ 



2005 



Fortney & Marleyl l2007h . different models con- 



verge toward the same effect on the planet atmosphere 
and evolutionary properties. Atmospheric thermal pro- 
files are strongly modified by irradiati on effects (see Bur 
rows chapter, this vo lume, as well as iBarman et al. I l200l[ 



Sudarskv et al.ll2003b as illustrated in Fig. |9] (upper panel). 



The heating of the outer layers by the incident stellar flux 
yields to an isothermal layer between the top of the convec- 
tive zone and the region where the stellar flux is absorbed. 
The top of the convective zone is displaced toward larger 
depths, compared to the non-irradiated case. The main ef- 
fect of the shallower atmospheric pressure-temperature pro- 
file is to drastically reduce the heat loss from the planet's 



11 



interi or, which can main tain higher entropy for a longer 
time ( Guillot et alJ 19961) . Consequently, the gravitational 
contraction of an irradiated planet is slowed down com- 
pared to the non-irradiated counterpart and the upshot is a 
larger radius at a given age (see Fig. |9] lower panel) 

The quantitative effect on the planet's radius depends 
on the planetary mass, parent star properties and orbital 
distance. Typical effects on the radius of irradiated Sat- 
urn mass or Jupiter mass planet located at orbital distances 
ranging between 0.02 AU and 0.05 AU around a solar- 
type star are of the order of 10%-20% jBaraffe et al.ll2003l: 



iBurrows et al ?2003UChab rier et al.ll2004l: lArras & BildstenI 
i2006: Fortney et al. 2007). A consistent comparison be- 
tween the theoretical radius and the observed transit ra- 
dius requires an additiona l effect du e to the thickness of 



the p l anet a tmosphere (Ba raffe et al. 2003; B urrows et al 



20031 l2007h . The measured radius is a transit radius at a 



given wavelength, usually in the optical, where the slant 
optical depth reaches ^1. This is at atmospheric layers 



above the photosphere dHubbard et al.ll200ll iBurrows et al 



20031) . The latter region is defined by an averaged normal 
optical depth r ^1, where the bulk of the flux is emitted 
outward and which corresponds to the location of the the- 
oretical radius. The atmospheric extension due to the heat- 
ing of the incident stellar flux can be significant, yielding 
a measured radius larger than the simple theoretical radius. 
This e ffect can add a few % (u p to 10%) to th e measured 
radius JBaraffe et al.l E)03: Burr ows et al]|2003h . Effects of 
irradiation on both the thermal atmosphere profile and the 
measured radius must be included for a detailed compari- 
son with observations and can explain some of the less in- 
flated exoplanets. They are however insufficient to explain 
the largest radii of currently known transiti ng exoplanets, 
such as TrES-4 with a radius R — 1.78i?j JSozzetti et al 



20091) . This fact points to other mechanisms to inflate close- 



in planets. 

4.3 Determining Transiting Planet Composition 

The fracti on of planets that are larger than can easily 
be explained (iMiller et al. I l2009[ put his fraction at ~40%) 
serve as a handicap to the goal to understand the compo- 
sition of these planets. We would like to be able to con- 
strain the fraction of planetary mass that is the non-H/He 
heavy elements, and understand how this varies with plane- 
tary mass, stellar mass, stellar metallicity, distance from the 
parent star, and other factors. Work on understanding the 
heavy element enrichment of planets has progressed, but 
must be regarded as incomplete, due to the poorly under- 
stood large-radius planets. 

It is certainly clear that a number of exoplanets, like 
our solar system giant planets (see are enriched in 
heavy material, as these planets have radii smaller than pure 
H/He objects. This idea is supported by our current un- 
derstanding of planet formation via the core-accretion sce- 
nario. This model is itself consistent with current obser- 
vations showing that metal rich environments characterized 




-0.5 
log M/Mj„p 

Fig. 8. — Mass-radius diagram for the known transiting 
planets (data are taken from the web site of Frederic Pont: 
www.inscience.ch/transits). Three iso-density curves are shown 
for the mean densities of Saturn (p = 0.62), Jupiter (p = 1.24) 
and Neptune (p = 1.61). These three solar system giant planets 
are also indicated by solid points. 
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Fig. 9. — Effect of irradiation for a planet at 0.05 AU from the 
Sun. Upper panel: T — P profiles of atmosphere models with 
intrinsic Teff=1000 K (a high value representative of giant planets 
at young ages) and surface gravity log 3=3 (cgs). The solid line 
corresponds to the irradiated model and the dashed line to the non- 
irradiated model. The locations of the photosphere (r ^ 1) and of 
the top of the convective zone are indicated by symbols (Models 
after Barman et al., 2001). Lower panel: evolution of the radius 
with time of a 1 Mj giant planet. Irradiated case: solid line; Non 
irradiated case: dashed line (Models aftci Baraffe et al, 2003). 
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by the high metallicity of the parent star favors planet for- 
mation ( Udrv & Santosll2007h . 

The exoplanets with the largest fraction of mass in heavy 
elements known at the time this chapter is written are 
WASP-7b, \yith M „ = 0.96Mj and Rp = 0.915i?j 
dHellier et al l '2009V HAT-P-3b with Mp = 0.6A/j and 



Rp = OMRj (Torres et al. 2007) and HP 149026 b with 
Mp = 0.36Mj and Rp = 0.755i?j dSato et al.l2005h . Their 
small radii indicate a global mass fraction of heavy ma- 
terial greater than what Jupiter contains (more than 12%) 
and in the case of HD 149026b, significantly greater than 
what Saturn contains (more than 30%). This suggests that 
the presence of a significant amount of heavy material is a 
property shared by all planets. Evolutionary models must 
thus take into account such enrichment and many efforts 
are currently devoted to the construction of a wide range of 
models with different amounts of heavy elements. Explo- 
ration of the effects of materials of different composition is 
limited by the available equations of state under the con- 
ditions of temperature and pressure characteristic of giant 
planets. As a simple starting point, and given the large un- 
certainty on the nature of heavy elements and their distri- 
bution inside planets, current models often assume that all 
heavy elements are located in the core and are water/ice, 
rock and/or iron. A possible interior density profile of HD 
149026b with a core of either ice or rock is shown in Figure 
[TOl Uncertainties in current planetary models due to uncer- 
tainties in the available EOS, the distribution of heavy ele- 
ments, and their chemical compo sition have been analyzed 
in detail in i Baraffe et al.l (120081) . Depending on the total 
amount of heavy material and its composition, the radius of 
an enriched planet can be significantly smaller compared to 
that of H/He dominated giant planets. This is illustrated in 
Fig. [TTl which compares the radius evolution with time of 
planets with different core sizes (in M^) and different core 
compositions. 

As noted, the problematic planets are those that are "too 
large," such that we cannot constrain the heavy element 
mass in their interiors. However, if the "true" additional 
energy source or contraction-stalling mechanism affects 
nearly all hot Jupiters to some degree (which seems likely 
given that 1/3 to 1/2 of planets are too large), then we may 
be able to take out its effect on planetary radii, and examine 
the heavy element enrichments of the transiting planet s as a 
collection. This was first done by lGuillot et al. who 



postulated an unspecified additional interior energy source 
equal to 0.5% of the absorbed incident energy. All plan- 
ets then required at least some amount of heavy elements 
to match the measured radii, and they found a correla- 
tion between stellar m etallicity and planetar y heavy element 
mass. Subsequently, Burrows et al.l (120071) found a similar 



correlation, using only enhanced atmospheric opacities for 
all planets, although t his cannot exp lain the largest-radius 
planets. Importantly, IGuillotl (l2008h find that the stellar- 
metallicity / planetary heavy-element correlation continues 
to hold with a larger number of planets. This appears to 
fit well with expectations for the core-accretion theory of 



HD 149026b (rock core) 
HD 149026b (ice core) 
Neptune 
Saturn 
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Fig. 10. — Interior density as a function of normalized radius for 
two possible models for HD 149026b compared with Neptune and 
Saturn. All planet models have been norms Jized to the radius a t 
which P=l bar. The Neptun e profile is fromlPodolak et al] ( 1 19951) 
and the Saturn profile is from lGuillollh999l) . The Saturn and Nep- 
tune models have a two-layer core of ice overlying rock, but this is 
only an assumption. For Neptune in particular the interior density 
profile is uncertain. The two profiles of HD 149026b assume a 
metallicity of 3 times solar in the H/He envelope and a core made 
entirely of either ice or rock. (After Fortney et al. 2006.) 
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Fig. 11. — Effect of heavy element enrichment on the evolution 
of the radius with time of a 1 Mj giant planet. Solid line: H/He 
planet; dashed line: Mcorc = 63 (ice); long dashed line: 
Mcoro = 159 Me (ice); dotted line: Mcore = 159 Me (rock); 
dash-dotted line: Mcorc = 159 Me (iron). 
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giant planet formation dDodson-Robinson & Bodenheimei 
2009h . Determining whether this and perhaps other rela- 
tions hold in the future will be an important area of transit- 
ing planet research in the future, as the number of detected 
planets continues to climb. 

4.4 Proposed Solutions To The Radius Anomaly 

As mentioned in §3.1, a significant fraction of transiting 
planets have large radii that cannot be explained by the ef- 
fect of thermal irradiation from the parent star alone. Other 
mechanisms are required to explain these inflated planets. 
Whether these mechanisms are peculiar, operating only un- 
der specific conditions, for certain planets, or whether basic 
physics is missing in the modeling of close-in giant planet 
structure are still open questions. In either scenario the in- 
flation mechanism(s) must be common. As noted, the lat- 
ter alternative stems from the suspicion that all observed 
transiting planets may contain a certain amount of heavy 
material. The denser planetary structure resulting from the 
presence of heavy material could counteract the effect of 
the "m issing" mechanism, yielding an observed "normal" 
radius ( Fortney et al ]|2006). 

Many mechanisms have been proposed since the discov- 
ery of HD 209458b, but no mechanism has gained the con- 
sensus of the community as being clearly important. We 
describe below the main ideas and comment on their status. 

• Atmospheric Circulation 

Based on numerical simulations of atmospheric cir- 
culation on hot Jupiters, Showman and Guillot (2002) 
suggested a heating mechanism in the deep interior 
driven by strong winds blowing on the planetary sur- 
face. The idea is that stellar irradiation produces 
strong day-night temperature contrasts, which drive 
fast winds. Their simulations produced a downward 
kinetic energy flux, about 1% of the absorbed stel- 
lar incident flux, which dissipates in the deep layers. 
This mechanism provides an extra source of energy 
which slows down the planet evolution and can ex- 
plain the large observed radii {Guillot and Showman, 
2002; Chabrier et al. 2004). 

This attractive and original scenario still requires 
confirmation. The very existence of the downward 
kinetic energy flux and its strength strongly depend 
on the outcome of atmospheric circulation models. 
Although this field has almost half a century history 
for solar system planets, it is still in its infancy when 
applied to hot Jupiters where conditions are very dif- 
ferent. Current models show important divergences 
(see the chapter by Showman et al.) and their robust- 
ness to describe the conversion process of heat from 
the parent star into me chanical energy (i.e winds) has 
even been questioned ( Goodmanl2009 ). The substan- 
tia l downward transport of k inetic energy reported 
bv lShowman & GuillotI ( 120021) has not been found in 
somewhat similar numerical simulations performed 



bv lBurkert et al.l ( l2005h . Moreover, the physical pro- 
cess which converts kinetic energy into thermal en- 
ergy in the Showman and Guillot scenario still needs 
to be identified. 

Enhanced Atmospheric Opacities 



Burrows et al.l (l2007n suggest that abnormally large 



radii of transiting planets can be explained by invok- 
ing enhanced opacities. Enhanced opacities could for 
instance be due to missing or underestimated opaci- 
ties in the c urrent generation of m odel atmospheres. 
In practice, iBurrows et al.l (l2007h implement these 
higher opacities by simply using supersolar metallic- 
ities (up to 10 times solar). The main effect of en- 
hanced opacities is to retain internal heat and to slow 
down the contraction. These authors conclude that a 
combination of enhanced opacities and the presence 
of dense cores can provide a general explanation to 
the observed spread in radii of transiting planets. In- 
terestingly, this work also finds the same correlation 
as suggested by iGuillot et al. I (l2006h between planet 
core mass and stellar metallicity. 

If opacity calculations are indeed correct, then higher 
opacities could be due only to increased metallicity. 
This would, however, present several weaknesses. 
It may be difficult to maintain a substantial frac- 
tion of heavy elements in strongly irradiated atmo- 
spheres, which are radiative down to deep levels. 
Perhaps mixing via atmospheric dy namics suffices 



to keep the atmosphere well mixed (iShowman et al 



a.; Sp i egeLe t ak 2009). Moreover, as discussed 
in lOuiUorfeoOSl) . iBurrows et alJ (l2007l) do not take 



into account the subsequent increase of molecular 
weight due to the increase of the planet atmosphere 
metallicity. This effect counters the effect of en- 
hanced opacities and may even dominate in some 
cases, yielding the opposite effect on the radius (see 
discussion in Guillot, 2008). Finally, for the most in- 
flated transiting planets yet detected, models with 
enhanced opacities are unable to reproduce their 
radii and anoth er additional mechanism is required 
dLiu et alJl2008h . 



Tidal Effects 

Another alternative, first suggested bv lBodenheimer et al 
laooi), is that tidal forces may heat the planet. There 
has come to be a large body of work in this area. If 
a planet has an eccentric orbit or a non synchronous 
rotation, internal tidal dissipation within the planet 
produces an energy source, which can slow down its 
contraction, or re-inflate th e planet after a previous 
contraction phase (see, e.g. lGu et al.ll2003l) . Follow- 
ing this idea, Liu et al. (2008) suggest a combina- 
tion of enhanced atmospheric opacities, the presence 
of heavy elements and heating due to orbital tidal 
dissipation, assuming non-zero eccentricities, to ex- 
plain the most inflated planets. Almost all studies on 
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tidal effects have assumed a tidal equilibrium state, 
yielding very short timescales for synchronization 
(- 10'^ - 10^ yr) and circularization (- 10^ - 10^ yr) 
compared to the estimated age of known exoplane- 
tary systems. This assumption implies that for eccen- 
tricity to be maintained on several Gyr and contribute 
to the heating, it should be continually excit ed. In 



the case of HD 209458b. Bodenheimer etaP (l200lh 
suggested the presence of an unseen planetary com- 
panion, which could force eccentricity. 

Tidal heating due to a finite current eccentricity may 
explain some of the large transit radii currently ob- 
served. It is, however, certainly not the mechanism 
which could explain all inflated planets. Constraints 
on the eccentricity of HD 2094 58b based on the tim- 



ing of the secondary eclipse (iDeming et al.l 120051) 



yield that a non-zero eccentricity is very unlikely. 
This explanation seems also improbab le for TrES- 
4 based on Spitzer observations by Knutson et al 



(l2009h who can rule out tidal heating at the level re- 
quired by Liu et al. (2008) to explain this planet's 
bloated size. 



More recentlv. lJackson et al.l(l2008h and lLevrard et al. 
(i2009l) have revisited the tidal stability of exoplanets. 
Essentially all examined transiting planets have not 
reached a tidal equilibrium state, implyin g that they 



wiU u lti mately fall onto the ce ntral star dLevrard et al 



2009h . iJackson et alj I (l2009h . confirming this view 
find that it is the youngest parent stars that tend to 
harbor the closest-in hot Jupiters, implying a loss of 
close-in planets with time. More importantly, these 
works stress that conventional circularization and 
synchronization timescales, which are widely used in 
the community, are i n most cases not correct. In the 
Levrard et al.l (l2009l) view, nearly circular orbits of 
planetary transiting systems currently observed may 
not be due to tidal dissipation. Conversely, obser- 
vations of non-zero eccentricity would be naturally 
explained without the need for gravitational interac- 
tions from undetected companions. 

Ibgui & Burrowsl ( l2009h and lMiller et all ( l2009l) have 



coup led a standard 2nd- order tidal evolution theory 
(e.g. Jacks on et al.ll2009h to planet structural evolu- 



tion models for close-in giant planets in single-planet 
systems to investigate under what circumstances tidal 
heating by recent eccentricity damping (which leads 
to an energy surge and radius inflation), together with 
tida l semimaior axis de cay, could inflate these plan- 



Miller et all (12009') find that this radius inflation 



ets. 

can occur for some planets at Gyr ages, perhaps ex- 
plaining some large radii, but this mechanism is un- 
likely to explain all inflated radii. 

Tidal heating could also be produced by a large obliq- 
uity, i.e the angle between the plane tary spin axis and 
the orb ital normal, as suggested bv IWinn & Hoknan 
(120051). Since planet obliquity is expected to be 



rapidly damped by tidal dissipation, a possibility to 
maintain a non-zero obliquity is to be locked in a 
Cassini state, i.e a resonance between spin and or- 
bital precession. However Levrard et al. (2007) show 
that although the probability of capture in a spin- 
orbit resonance is rather good around 0.5 AU, it de- 
creases dramat i cally with semi-major axis. Also, 
iFabrycky et alJ ( 2007 ) rule out possible drivers, such 
as the presence of a second planet, of a high obliquity 
Cassini state for HD 209458b. They conclude that 
very special configurations are required for obliquity 
tides to be an important source of heating. On the 
whole, it thus seems difficult to invoke tidal dissipa- 
tion as the mechanism which could explain all cur- 
rently observed inflated exoplanets, although impor- 
tant quantitative effects have been identified. Consid- 
erable work in this area continues. 

Double Diffusive Convection 

Chabrier & Baraffe (120071) uggest that the onset of 



layered or oscillatory convection, due to the presence 
of molecular weight gradients, can reduce heat trans- 
port in planetary interiors and slows down the con- 
traction, providing an explanation for the large spread 
in radii of transiting planets. The formation of lay- 
ers is a characteristic of double-diffusive convection 
which may occur in a medium where two substances 
diffuse at different rates. This is a well known process 
in oceans or salty lakes where the two substances are 
heat and salt. Applied to the interior of planets, in the 
presence of a compositional gradient, convection can 
break into convective layers separated by thin diffu- 
sive layers. The heat transport efficiency is thus sig- 
nificantly reduced, because of the presence of multi- 
ple diffusive layers, compared to the case of a fully 
homogeneous planet where convection is assumed 
to be fully adiabatic. This proces s is similar to th e 
so-called semiconvection in stars (^ Stevensoiil 19791) . 
Based on a phenomenological approach, and assum- 
ing a molecular weight gradient in the most inner part 
of the planet, Chabrier and Baraffe (2007) show that 
a significant number of diffusive layers strongly re- 
duce the heat escape. They suggest that the composi- 
tion gradient is inherited from the formation process, 
during accretion of planetesimals and gas, or due to 
core erosion. The upshot is a significantly inflated 
planet compared to its homogeneous and adiabatic 
counterpart. 

The idea that planetary interiors may not be com- 
pletely homogeneous and convection not fully ef- 
ficient was already suggested by D. Stevenson for 
Jupit er |Stevenson 1985) and for Uranus and Nep- 



tune (iHubbardet al.lll995l) . (See §3.) Double diffu- 



sive convection is indeed a well known process under 
Earth conditions, which have similarities with those 
found in the interior of giant planets (see Chabrier 
and Baraffe, 2007 for details). Whether this scenario 
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can explain all inflated transiting planets is still de- 
batable. The key questions are whether the diffusive 
layers commonly form, and can survive on timescales 
of Gyr, characteristic of the age of the transiting ex- 
oplanets. Development of 3D numerical simulations 
of this process under planetary conditions could pro- 
vide clues on its long term existence in giant planet 
interiors. 

4.5 Hot Neptune Planets 

The dream of discovering extrasolar planets of a few 
Earth masses is now becoming reality. More than a 
dozen planets with masses in the Uranus-Neptune range 
( " <: 20 Mffi) h as been detected by radial velocity surveys 
(|Udry & San tos 2007). Because these light planets are very 
close to the detection threshold, their discovery suggest 
that they are rather common. Given their low mass and 
close orbit, the question of their origin deserves some at- 
tention. Formation models based on the core accretion 
scenario suggest that hot Neptunes are composed of a 
large heavy material core (ice or rock) and have formed 
without accumulation of a substantial ga seous envelope 
jAlibert et al.ll2006t iMordasini et al.ll2008l) . Another sug- 
gestion is that hot Neptunes could have formed from more 
massive p rogenitors and hav e lost most of their gaseous 
envelope te araffe et al.ll2005h . The latter idea arises from 
observational evidences that close-in exoplanets may un- 
dergo evaporation proc esses induced by the hig h energy 
flux of the parent star (iVidal-Madjar et al.ll2003h . as well 
as high mas s- loss r ates found in some early models (e.g. 



Lammer et alfcoosl) 



The interpretation of the IVidal-Madiar et al.l (l2003h ob- 
servations of an extended neutral hydrogen atmosphere 
around HD 209458b is an active area that is still controver- 



greater than 85%, which is the overall heavy element con- 
tent of Uranus and Neptune (Guillot, 2005). Given current 
uncertainties on planetary interior structure models, only 
the bulk of heavy elements can be inferred for each, about 
20 Me for a total mass of 22 M® for GJ 436b, which 
is the more well-studied of the two (see Fig. [T2] |. Although 
small in terms of mass, the contribution of the H/He en- 
velope to the total planetary radius is significant. This is 
illustrated in Fig. [12] where the radii of pure water and pure 
rocky planets of the same mass are also indicated. 

Better determination of the chemical composition of 
heavy material and its distribution within the planet must 
await improved EOSs for H/He and heavy elements 
(Baraffe et al, 2008). These discoveries however confirms 
the large heavy element content that can be expected in ex- 
trasolar planets and supports the general picture of planet 
formation drawn by the core accretion model. 

4.6 Young Giant Planets 

Giant planet thermal evolution models are being tested 
at Gyr ages for solar system planets and the transiting plan- 
ets. It is clear from giant planet formation theories (see the 
chapter by D' Angelo & Lissauer) that these planets are hot, 
luminous, and have larger radii at young ages, and they con- 
tract and cool inexorably as they age. However, since the 
planet formation process is not well understood in detail, 
we understand very little about the initial conditions for the 
planets' subsequent cooling. Since the Kelvin-Helmholtz 
time is very short at young ages (when the luminosity is 
high and radius is large) it is expected that giant planets 
forget their initial conditions quickly. This idea was estab- 
lished with the initial Jupiter cooling models in the 1970s 



(IGraboskeet al.l ll975: Bodenheimer 1976). 



sial ('Ben-Jaffel'2007^.' Holmstrom et alJ2008l:IVidal-Madjar et 
I2OO8; Ben-Jaffel 20081). However, at flie same time, the 
various groups computing models of evaporative mass loss 
at small orbital distances have been converging to mass 
loss rates that yield a total mass loss of only ^1% for 
HP 209458b ov e r die lifetime of die s ystem (Yelle 2004; 
iTian et al.1 12005"; 'Gar cia Mufio^ 120071: Ivelle et al., ,2008,: 
Murrav- Clay et al. 200^ If this is correct, than these evap- 
oration processes are too small to significantly affect the 
evolution of the currently observed close-in planets. Small 
evaporation rates seem also to be c onsistent with the ob - 
served mass function of exoplanets jHubbard et al. 2007). 
Improved statistics in the low planetary mass regime, con- 
firmation of the observations of Vidal-Madjar et al. (2003), 
and the extension of similar observations to other transiting 
planets will shed light on these issues. 

Independently of these issues, the interior properties of 
two extrasolar hot Neptune planets were recently revealed 
by the remarkable disco very of the fi rst transiting Neptune- 
mass planets GJ 43 6b dGillon et alj r2007) and HAT-P-llb 
dBakos et al. I l2009bh . Both planets have a radius compara- 
ble to that of Neptune, indicating heavy material enrichment 



Since our solar system's giant planets are thought to be 
ail 5 Gyr old, there is little worry about how thermal evolu- 
tion models of these planets are affected by the unknown 
initial conditions. The same may not be true for very young 
planets, however. Since giant planets are considerably 
brighter at young ages, searches to directly image planets 

these searches are 



now bearing fruit fchauvin et al. 


2OO5I 


Marois et al. 2008; 


Kalas et al. 2008; Laaranae et al. 


20091). It is at ages of a 



few million years where understanding the initial conditions 
and early evolution history is particularly important. Tra- 
ditional evolution models (|2J, which are applied to both 
giant planets and brown dwarfs, employ an arbitrary start- 
ing point. The initial model is large in radius, luminosity, 
and usually fully adiabatic. The exact choice of the starting 
model is often thought to be unimportant, if one is inter- 
ested in following the evolution for ages greater than 1 Myr 
dBurrows et al.lll997l:lchabrier & Baraffe 2000). 

Thermal evolution models, when coupled to a grid of 
model atmospheres, aim to predict the luminosity, radius, 
Tcff , thermal emission spectrum, and reflected spectrum, as 
a function of time. When a planetary candidate is imaged, 
often only the apparent magnitude in a few infrared bands 
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are known. If the age of the parent star can be estimated (it- 
self a tricky task) then the observed infrared magnitudes can 
be compared with calculations of model planets for various 
masses, to estimate the planet's mass. Recall that mass is 
not an observable quantity unless some dynamical informa- 
tion is also known. It is not known if these thermal evolu- 
tion models are accurate at young ag es-they are re l atively 
untested, which has been stressed by Ba raffe et alJ (l2002h 
fo r brown dw arfs and' Marley et al. ( 2007h for planets. 



Marie v et al.. (J007) examined the issue of the accuracy 



of the arbitrary initial conditions (termed a "hot start" by the 
authors) by using initial conditions for cooling that were 
not arbitrary, but rather were g iven by a leading co re ac- 
cretion planet formation model ( Hubickyj et aljEo05 ). The 
core accretion calculation predicts the planetary structure at 
the end of f ormation, wh e n the planet has reached its final 
mass. The iMarlev et al I (l2007b cooling models used this 
initial model for time zero, and subsequent cooling was fol- 
lowed as in previously published models. Figure [13] shows 
the resulting evolution. The cooling curves are dramatically 
different, yielding cooler (and smaller) planets. The ini- 
tial conditions are not quickly "forgotten," meaning that the 
cooling curves do not overlap with the arbitrary start models 
for 10^ to 10^ years. What this would mean, in principle, 
is that a mass derived from "hot start" evolutionary tracks 
would significantly underestimate the true mass of a planet 
formed by core accretion. 

Certainly one must remember that a host of assump- 
tions go into the formation model, so it is unlikely that 
these new cooling models are quantitatively correct. How- 
ever, they highlight that much additional work is needed 
to understa nd the energeti c s of t he planet formation pro- 
cess. The iHubickv) et al. I (l2005 l) models yield relatively 
cold initial models because of an assumption that accret- 
ing gas is shocked and readily radiates away this energy 
during formation. This energy loss directly leads to a low 
luminosity starting point for subsequent evolution. Signif- 
icant additional work on multi-dimensional accretion must 
be done, as well as on radiative transfer during the accretion 
phase, before we can confidently model the early evolution. 
Thankfully, it appears that detections of young planets are 
now beginning to progress quickly, which will help to con- 
strain these models. 



5. FUTURE PROSPECTS 



The future of understanding the structure, composition, 
and evolution of giant planets is quite promising. Most im- 
mediately, and with the biggest impact, will be the detec- 
tion of more Neptune-class transiting planets, in addition 
to GJ 436b and HAT-P-llb. Of particular interest will be 
the mass-radius relation of planets around ~10 AI(^, since 
this is estimated to be a boundary between planets that have 
H/He envelopes and those that do not. A radius that is larger 
than that calculated for a pure water planet is unambiguous 




Fig. 12. — Evolution of a planet with the mass of GJ 436b (22.6 
M©) and different heavy element contents and compositions. The 
solid line corresponds to a model with a water core of 21 M© and 
the dashed line to a rocky core of 19.5 (after the models of 
Barqjfe et al, 2008). The long dashed line indicates the radius of 
a pure water planet and the dash-dotted line corresponds to a pure 
rocky planet (after the models of Fortney et al. , 2007) 




Fig. 13. — T hermal evolut i on of giant planets from 1 to 10 Mj, 
adapted from iMarlev et al.l ( l2007b . The dotted curves are stan- 
dard "hot start" models with an arbitrary initial condition, and the 
solid curv es use as an initial cond ition the core accretion formation 
models of lHubickvi et alj ( l2005h . 
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Fig. 14. — Contraction of a 1 Mj planet, with 25 Af® of heavy 
elements (in a core), over time. The planets are placed at 0.02, 
0.045, 0.1, 1.0, and 10 AU from a constant luminosity Sun, to 
sho w the effect of s t ellar i rradiation on thermal evolution. Adapted 
from lFortnev et alj j2007h . 



evidence of a H/He envelope. Certainly the diversity of radii 
for the transiting gas giants has been surprising, and under- 
standing Neptune-class radii as a function of mass and or- 
bital distance will be fascinating. 

Several ground-based surveys for transiting planets are 
now scouring the northern hemisphere for Jupiter-class 
planets, with some additional attention being paid to the 
southern hemisphere. The required photometric precision 
to detect a Jupiter-type transit (a 1 % dip in stellar flux) is not 
difficult to achieve. Furthermore, the detection of HAT-P- 
1 lb opens up the possibility that these same ground-based 
surveys will also be able to detect significant numbers of 
Neptune-type planets as well. In addition, the implemen- 
tation of orthogonal transfer array CCDs may allow for 
ground-based photometric pr ecision that ap proaches that 
of space-based platforms (e.g. Ijohnson et al.ll2009.) . which 
would significantly help the cause of detecting smaller plan- 
ets from the ground. 

The CoRoT mission, which has already announced five 
planets, and will continue to search for a few more years, 
will be important in adding to our sample size in two ways. 
The first is of course smaller planets. The second is increas- 
ing the sample size of transiting planets at wider orbital sep- 
arations. CoRoT is surveying particular areas of the sky for 
120 days, meaning that 40-day orbits of giant planets can 
be easily confirmed via 3 detected transits. Around a solar- 
type star, a 40-day orbit yields a semi-major axis of 0.23 
AU, which is a factor 33 reduction in flux from that inter- 
cepted at 0.04 AU. Understanding radii at the largest possi- 
ble range of incident fluxes may allow us to understand the 
reason(s) for the large radii of the currently detected plan- 
ets. 

The Kepler mission, which began taking science data in 
May 2009, is an even more ambitious telescope to detect 
smaller transiting planets in longer period orbits. Although 



its main goal is to ascertain the frequency of Earth-radius 
planets in Earth-like orbits around Sun-like stars, it will also 
do the same for larger planets. Kepler will be able to detect 
multiple transits of planets out to 1 AU, a factor of over 600 
reduction in flux from 0.04 AU. In addition, even longer 
period transiting giant planets may be detected, if suitable 
followup is done. Since some mechanisms proposed to ex- 
plain the large radii of the close-in planets should be signif- 
icantly muted at larger orbital separations, finding planets 
farther from their parent stars will likely be the most impor- 
tant step in understanding what leads to these large radii. In 
Figure [14] we show a specific prediction for the contraction 
of a 1 A/j planet over a factor of 250,000 i n incident stel- 
lar flux. As detailed in lFortney et alJ (l2007h . the effects of 
stellar flux are muted beyond the current group of close-in 
transiting planets, and planets out to ^1 AU should have 
radii quite similar to those at ~0. 1 AU, in the absence of 
missing physics for these more distant planets. 

The orbital dynamics of particular exoplanets in some 
systems may give us direct constraints on a planet's inte- 
rior state. The apsidal precession rate of a planetary or- 
bit is directly proportional to the tidal Love number, /c2. 
This number parameterizes the internal density distribu- 
tion of a planet. Authors have recently pointed out in- 
stances in which k2 could be measured or well-constrained. 
These include precession due to a tidal-induced gravita- 
tional quadrupole on the planet by its parent star, for very 
close-in planets (Ragozzine & Wolf 2009), which could be 
measured as a change in transit shape over time. Another af- 
fects transiting planets in multi-planet systems, which now 
only includes HAT-P-13b (Bakos et al. 2009a). For this sys- 
tem a refined measurement of current planetary eccentricity 
can constrai n k-2 as w ell as Q, its tidal dissipation quality 
factor (Baty gin et aDE )09). 

Detections of transiting giant planets at younger ages 
would be very important as these planets would inform our 
understanding of contraction with time in the face of intense 
stellar iiTadiation. This would shed light on the initial con- 
ditions for evolution, post-formation, as well as allow us to 
better understand the nature of the physical process that is 
causing large radii at gigayear ages. Towards this goal, sev- 
eral transit surveys of open clusters have been performed, 
but they have not netted any planets to date. 

The focus of comparing models to observations is al- 
ready shifting from specific planets to samples of planets, 
and will soon shift to a statistically significant number of 
planets, with the additional detections from the ground and 
from space. A good reference for t he kin d of work, just 
starting to be done, is iFressin et al. ( 2007 ). who analyzed 
in detail the OGLE transiting planet survey. They simulate 
the OGLE survey: given the properties of the thousands of 
stars that were monitored, the known planet frequency as a 
function of stellar mass and of orbital distance, and imple- 
ment giant planet contraction models, to derive constraints 
on, for instance, possible separate populations of plane ts. 

The planets dire ctly imaged by iMarois et al.l (120081) and 
Kalas et al.l (120081) have fully opened the door to direct 
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imaging, whic h began yielding plan etary candidates a few 
years ago (e.g. IChauyin et al.ll2005h . The characterization 
of these planets will present different challenges compared 
to the transiting planets. While transiting planets yield ac- 
curate masses and radii, the atmospheric characterization is 
challenging. For directly imaged planets, masses and radii 
likely cannot be measured, but spectra should be more eas- 
ily obtained. Spectra can yield the planet's Tcff, which can 
be directly compared to thermal evolution models. How- 
ever, current techniques are limited to planet-to-star flux 
ratios of ^ 10~^. New instruments coming online in the 
very near future, such as the Gemini Planet Imager (GPI) 
on Gemini South and Spectro-Polarimetric High-contrast 
Exoplanet REsearch (SPHERE) at the VLT will allow for 
contrasts of ~ 10~^ — 10^^. This will allow for the direct 
imaging and characterization of giant planets at a variety of 
masses, and also a variety of ages, so the first ~ 1 — 100 
Myr of giant planet evolution should be relatively well un- 
derstood via observations. 

We are now over a decade into the new era of studying 
giant planets as a class of astronomical objects. It is the 
expansion of this class beyond Jupiter, Saturn, Uranus, and 
Neptune, that will enable us to better understand the forma- 
tion and evolution of these planets. Much like our under- 
standing of the stars is greatly enhanced by studying more 
than just the Sun, so will our understanding of giant planets 
grow. 

Over the past several years we have seen strange and 
startling transiting planets with hugely inflated radii, and 
those with small radii that must include vast interior stores 
of heavy elements. They have expanded our imaginations 
regarding the possible structure of giant planets. We are 
now beginning to see these original oddballs within the con- 
tinuum of giants planets that extend far beyond what we see 
in our solar system. Astronomers will do doubt continue to 
creatively find ways to detect and characterize many more 
planets in the future. 
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